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This Letter describes a model-independent search for production of new resonant states in pho- 
ton+jet events in 2.11 fb _1 of proton-proton collisions at ^/s — 7 TeV. We compare the photon+jet 
mass distribution to a background model derived from data and find consistency with the background 
hypothesis. Given the lack of evidence for a signal, we set 95% CL limits on generic Gaussian-shaped 
signals and on a benchmark excited-quark (q*) model, excluding 2 TeV Gaussian resonances with 
cross section times branching fraction times acceptance times efficiency near 5 fb and excluding q* 
masses below 2.46 TeV, respectively. 



PACS numbers: 12.60.-i,13.85.Rm,14.70.Bh,14.80-j 

In the Standard Model (SM), proton-proton (pp) col- 
lisions do not produce photon+jet pairs through a res- 
onance. Direct photon+jet production occurs at tree 
level via Compton scattering of a quark and a gluon 
or through quark-antiquark annihilation, the former ac- 
counting for the majority of direct photon+jet produc- 
tion at all center-of-mass energies. Events with a high 
transverse momentum photon and one or more jets can 
also arise from radiation off final-state quarks or from 
multi-jet processes where dijet or higher-order events 
produce secondary photons during fragmentation of the 
hard-scatter quarks and gluons [lj-|4j]. The photon+jet 
invariant mass (w 7J ) || distribution resulting from this 
mixture of processes is smooth and rapidly falling, con- 
stituting a promising place to look for resonances. 

Despite many possible exotic production mechanisms 
such as excited quarks H. 3, q uirks [81-Hoj]. Regge ex- 
citations of string theory [Ill4l4|. and topological pions 
[l5j |. the most recent searches for photon+jet resonances 
were published a decade ago [16l4l9l |. The previous most 
sensitive search for new phenomena in the photon+jet 
final state places limits on effective signal cross sections, 
cross section times branching fraction times acceptance 
times efficiency, of the order of 1 pb and on excited-quark 
masses up to 460 GeV at the 95% confidence level 161. 

This Letter describes a general search for resonant s- 
channel photon+jet production in 2.11 fb -1 of pp col- 
lisions at a center-of-mass energy y/s = 7 TeV with 
the ATLAS detector. It follows earlier measurements of 
isolated photon differential cross sections at the Large 
(LHC) 
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The entire m 7 j distri- 



Hadron Collider 
bution is fit to a smooth function to obtain the back- 
ground to this search. We look for evidence of a narrow 
resonance, not much wider than the detector mass resolu- 
tion. This search extends the method used in the search 



for resonant dijet production [22j to handle the more di 



verse mixture of processes contributing to the m 7J - distri- 
bution. In the absence of signal, we use Bayes' theorem 
to set limits on Gaussian-shaped resonances and on a 
benchmark excited-quark (q*) model @, 01- 



In the excited-quark model studied here, the LHC 
could produce single q* states with vector-like couplings 
to the W and Z gauge bosons via the absorption of a 
gluon by a quark. As in Ref. [zj, we define the model 
by one parameter, the excited-quark mass m q * , setting 
the compositeness scale equal to m 9 » and SU(3), SU(2), 
and U(l) coupling multipliers f s =/ = /' = 1. At 
Trig* = 2.5 TeV, this gives branching fractions for u* — > ug 
and u* — > 1/7 of 0.85 and 0.02, respectively. The corre- 
sponding branching fractions for d* quarks are 0.85 and 
0.005, respectively. We do not make any further assump- 
tions about higher-order corrections or the excited-quark 
dynamics and neglect scale uncertainties and uncertain- 
ties on parton distribution functions (PDFs) in order to 
provide a convenient benchmark process for theoretical 
reinterpretation. 

We simulate the SM direct photon processes and the 
q* model with Pythia 6.4.25 [H] using the AMBT1 
0, MSTW2007 PDFs [H, and a GEANT4-based 



tune 

detector simulation [26J, |27(. Supplementary studies of 
the background shape function are performed with next- 
to-leading-order JETPHOX 1.3.0 0,0. Additional inelas- 
tic pp interactions, termed pile-up, are included in the 
event simulation, distributed so as to reproduce the num- 
ber of collisions per bunch crossing in the data. The mean 
number of pile-up interactions is approximately 6. 

A detailed description of the detector is available in 
Ref. [lil]. Photons are detected by a lead/liquid- argon 
sampling electromagnetic calorimeter (EMC) with an ac- 
cordion geometry. In front of the EMC, the inner de- 
tector allows an accurate reconstruction of tracks from 
the primary pp collision point and also from secondary 
vertices, permitting an efficient reconstruction of photon 
conversions in the inner detector up to a radius of 80 cm. 
For | r\ | < 1.37 [29], an iron/scintillator-tile calorimeter 
behind the EMC provides hadronic coverage. The endcap 
and forward regions, 1.5 < \t)\ < 4.9, are instrumented 
with liquid-argon calorimeters for both electromagnetic 
and hadronic measurements. For further details relevant 
to photon identification and measurement, see Ref. [20j |. 
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For details relevant to jet detection and measurement, 
see Ref. Events are collected with a trigger requir- 
ing at least one photon candidate with transverse mo- 
mentum (pt) above 80 GeV. The i nteg rated luminosity 



32j,|33|- 



of the sample is (2.11 ± 0.08) fb" 1 

Events containing at least one photon and at least one 
jet are selected for analysis. Each event must have a 
primary vertex with at least five charged-particle tracks 
with pt > 400 MeV. Multiple vertices can appear when 
pile-up interactions occur for the same bunch crossing. 
If more than one vertex is found, the primary vertex is 
taken as the vertex with the highest scalar sum p\ of 
associated tracks. Photon candidates with pr > 85 GeV 
and \r\\ < 1.37 and jet candidates with px > 30 GeV and 
\r]\ < 2.8 are used. These objects are identified using cri- 
teria that closely follow those ap plie d in the isolated pho- 
ton cross section measurement [20j and dijet resonance 
search [22J. Subleading photons or jets are allowed; when 
more than one photon or jet is found, the highest pt can- 
didates are selected to constitute the photon+jet pair. 
Jets are reconstructed from clusters of calorimeter cells 



using the anti-/c t clustering algorithm [34[ with radius 
parameter R = 0.6. Jet energies are corrected to the 
hadronic scale [3(J 31 1. Jet candidates are rejected in 



regions of the calorimeter where the jet energy is not 
yet measured in an optimal way. Candidates consistent 
with spurious calorimeter noise or energy spikes are also 
rejected. 

Photon candidates are reconstructed from clusters in 
the electromagnetic calorimeter and tracking information 
provided by the inner detector. They satisfy standard 
ATLAS selection criteria that are desi gne d to reject in- 
strumental backgrounds from hadrons [20(. The photon 
candidates must meet pt- and 77-dependent requirements 
on hadronic leakage, shower shapes in the electromag- 
netic strip layer, and shower shapes in the second sam- 
pling layer of the electromagnetic calorimeter. Inner de- 
tector tracking information is used to reject electrons and 
to recover photons converted to e + e~ pairs. Energy cali- 
brations are applied to photon candidates to account for 
energy loss in front of the electromagnetic calorimeter 
and for both lateral and longitudinal leakage. Events 
are discarded if the leading photon appears in calorime- 
ter cells affected by noise bursts or transient hardware 
problems. 

These photon identification criteria reduce instrumen- 
tal backgrounds to a negligible level, but much of the 
substantial background from secondary (jet fragmenta- 
tion) photons remains. We reduce this background with 
requirements on nearby calorimeter activity. Associated 
"isolation" calorimeter energy near the photon candidate 
is calculated by summing the transverse momentum as 
measured in electromagnetic and hadronic calorimeter 
cells inside a cone of radius AR = y/ (Aij) 2 + (A</>) 2 = 
0.4 centered on the photon cluster, but excluding the en- 
ergy of the photon cluster itself. The isolation energy 



is corrected on an event-by-event basis for the ambient 
energy density due to pile-up and the underlying event. 
This isolation energy is required to be less than 7 GeV. 

The photon deposits energy in the electromagnetic 
calorimeter in such a way as also to be reconstructed as 
a jet. Jets within AR < 0.2 of the photon are therefore 
not considered in this analysis. We require an angular 
separation Ai?(7,jet) > 0.6 between the signal photon 
and all jets with pt > 30 GeV to reduce the background 
from photons during fragmentation of final state quarks 
(fragmentation photons) and to reduce the systematic ef- 
fects from leakage of nearby jet showers into the photon 
isolation energy measurement. 

Additional reduction of fragmentation photon back- 
ground is achieved by requirements on the photon and jet 
pseudorapidities. Dijet production rates increase with jet 
pseudorapidity whereas rates for our assumed s-channel 
signal would diminish. We restrict the analysis to pho- 
tons in the barrel calorimeter, |7/| < 1.37, and require 
1% — Vj\ < 1-4 between the photon and jet. The for- 
mer criterion was chosen to avoid kinematic bias of the 
?7i 7 j distribution due to inclusion of any r\ range where re- 
construction efficiency is lower, such as the barrel-endcap 
transition region 1.37 < |?7| < 1.52. The latter was chosen 
by optimizing expected significance using the |?7 7 — r\j \ dis- 
tributions found in excited-quark signal simulation and 
background-dominated control data selected as in the 
nominal analysis but inverting the photon isolation re- 
quirement. This control sample is also used to check the 
background estimate. 

After the above selections, Fig. Q] shows the distribu- 
tion of the m 7J - invariant mass in bins equal to the mass 
resolution. The m 7J - resolution is about 4% at 600 GeV, 
improving to 3% at 2 TeV. We determine the combined 
SM and instrumental background to the search by fitting 
this distribution to the four-parameter ansatz 
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The motivation for this function is discussed in Ref. IK 



35j-l37j. The fit result is also shown in Fig.[T] The bottom 
panel of the figure shows the statistical significance of 
the difference between data and the fit in each bin [38j |. 
With a negative log-likelihood test statistic, the p- value is 
23%, indicating the data distribution is compatible with 
Eq. [U The functional form also describes the leading- 
order Pythia direct photon prediction for comparable 
event statistics. 

We search for statistical evidence of a resonance in this 
distribution using the BumpHunter algorithm [39| . The 
algorithm operates on the binned m 7 j distribution, com- 
paring the background estimate with the data in mass 
intervals of varying contiguous bin multiplicities across 
the entire distribution. For each interval in the scan, it 
computes the significance of any excess found. The algo- 
rithm identifies the interval 784-1212 GeV, indicated by 
the vertical lines in Fig. [TJ as the single most discrepant 
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FIG. 1. Invariant mass of the photon+jet pair for events pass- 
ing the final selection. Overlaid: the fitted background func- 
tion integrated over each bin (stepped solid line), the most 
discrepant region identified by BumpHunter (two dashed 
vertical lines), and three examples of excited-quark signals, 
normalized to luminosity, as described in the text. The bot- 
tom panel shows the statistical significance of the difference 
between data and background in each bin. 



interval. The significance of the outcome is evaluated 
using the ensemble of possible outcomes for the signifi- 
cance of any region in the distribution in the background- 
only hypothesis, obtained by repeating the analysis on 
pseudodata drawn from the background function. Be- 
fore including systematic uncertainties, the probability 
(p- value) of observing a background fluctuation at least 
as significant as the above, including the trials factor, or 
"look-elsewhere" effect, is 20%. Inclusion of systematic 
uncertainties renders the p- value similarly large. 

Lacking evidence of any signal, we exclude two types 
of photon+jet signals: a generic signal with Gaussian 
distribution and arbitrary production cross section, and 
the excited-quark model. We compute Bayesian limits at 
95% credibility level (CL) using a prior probability den- 
sity that is constant for positive values of the signal pro- 
duction cross section and zero for unphysical, negative 
values, as described in Ref. (lb] . We consider system- 
atic uncertainties on expected signal yield due to imper- 
fect knowledge of the detector: the integrated luminosity 
(3.7%), trigger efficiencies (<0.5%), and signal photon 
identification efficiencies (2.0%). The last of these con- 
sists of isolation (0.4%), pile-up interactions (0.5%), con- 
versions (1.2%), simulation mismodeling (1.3%), and the 
extrapolation of the photon identification efficiency to 
high pt (<0.3%). Uncertainties on photon energy scale 



(0.5-1.5%), jet energy scale (2-4%), and jet energy res- 
olution (5-15%) contribute through their effects on the 
signal distribution. These systematic uncertainties are 
treated as marginalized Gaussian nuisance parameters in 
the limit calculation. 

We also evaluate two systematic uncertainties on the 
background estimate. To account for the statistical un- 
certainties on the background fit parameters, we repeat- 
edly fit the background function to pseudodata for each 
bin drawn from Poisson distributions. The mean of the 
Poisson distribution for a given bin corresponds to the 
number of entries actually observed in that bin in the 
data. We then take the variation in the fit predictions 
for a given bin, 0.5% of the background at low mass to 
almost 10% of the background at 2 TeV, as indicative 
of the systematic uncertainty. We treat this bin-by-bin 
uncertainty in the limit as fully correlated, using a sin- 
gle nuisance parameter that scales the entire background 
distribution. 

While our function can describe the m 7 j shape for di- 
rect photon production, as modeled in the Pythia di- 
rect photon+jet simulation, the function need not remain 
a good description of the full distribution after includ- 
ing nonisolated and fragmentation photon events. For 
example, the function describes the next-to-leading or- 
der prediction implemented in jetphox, which includes 
the fragmentation photon contributions, for some viable 
choices of theory parameters but not for others. 

The second background systematic uncertainty ac- 
counts for any unmodeled features of fragmentation pho- 
ton events in our isolated photon sample. We fit the 
background function to the m 7J - distribution in the con- 
trol data selected with the inverted isolation requirement, 
then measure for each m-yj bin the magnitude of any devi- 
ation from the fit, and assign the ratio of the deviation to 
the fit expectation as a parametrization bias systematic 
uncertainty. To extrapolate this uncertainty to large m^j 
where few control data exist, we fit the tail m 7J > 1 TeV 
with a two-degree polynomial. 

Figure [2] shows the model-independent limits on the ef- 
fective cross section, cross section a times branching frac- 
tion B times acceptance A times efficiency e, of a poten- 
tial signal as a function of the central mass of each signal 
template. We take the signal lineshape to be a Gaussian 
distribution with one of three widths, uqItiiq — 5%, 7%, 
and 10% of the central mass of the Gaussian. The limit 
weakens as the width increases and the peak becomes less 
distinct. For example, for a 1 TeV signal the limit for a 
width of 10% is 1.6 times the limit for a width of 5%. 

The limit on the effective cross section in the excited- 
quark model is shown in Fig. [3] as a function of the q* 
mass. Also shown are ±lcr and ±2cr uncertainty bands 
indicating the underlying distribution of possible limit 
outcomes in the background-only hypothesis. The solid 
line indicates the prediction from the Pythia excited- 
quark implementation. We exclude such excited quarks 
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with masses below 2.46 TeV at 95% CL, complement- 
ing the more stringent exclusion below 2.99 TeV on this 



specific q* model in the dijet final state [22 1 
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FIG. 3. The 95% CL upper limits on a X B x A x e for excited 
quarks decaying to a photon and a jet as a function of the 
signal mass, m q * , for the Pythia q* prediction and taking 
into account systematic uncertainties. 



FIG. 2. The 95% CL upper limits on a x B x A x e for a 
hypothetical signal with Gaussian m 7J distribution decaying 
to a photon and a jet as a function of the signal mass, mc, 
for three values of the relative Gaussian width ad tug and 
taking into account systematic uncertainties. 

In conclusion, the photon+jet mass distribution mea- 
sured in 2.11 fb^ 1 of pp collision data collected at y/s = 
7 TeV by the ATLAS Collaboration has been examined 
for narrow resonances. The observed distribution extends 
up to masses of about 2 TeV. It is well described by a 
smooth function fitted to it and assumed to represent the 
SM expectation. No evidence for the production of reso- 
nances is found. We set limits at 95% CL on Gaussian- 
lineshape and excited-quark signal using Bayesian statis- 
tics. The limits on Gaussian resonances, for example, ex- 
clude 2 TeV resonances with effective cross sections near 
5 fb. We also exclude the excited-quark model in the 
photon+jet final state for masses up to 2.46 TeV. The 
limits reported here on resonant production of new par- 
ticles in the photon+jet final state are the most stringent 
limits set to date in this channel. 
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